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Abstract It is shown that calculus can apply on a fractal structure with 
the condition that the infinitesimal limit of change of the variable is larger 
than the lower cut-off of the fractal structure, and an assumption called local 
decomposability. As an application, it is shown that the angular projection 
of a fractal distribution in 3-dimensional space is not homogeneous at suffi- 
ciently large angles. Therefore the angular projection of galaxy distribution 
for sufficiently large angles can discriminate the fractal and the homogeneity 
pictures. 
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It seems that calculus cannot apply in a fractal structure, which possesses 
infinite places of discontinuity. However, it is shown in this letter that it can 
be done in a conditional way. Just as that only with calculus can we deal 
with manifold and function with an arbitrary shape, the conditional calculus 
on fractals will enable us deepen the studies of fractal. 

This work is partly motivated by the investigations on fractal distribution 
of galaxies. While it is an agreement that the galaxy distribution approx- 
imates a fractal over a considerable range of scales, it is under a debate 
whether it becomes homogeneous on a scale about, say, 20h~^Mpc [0], 
or is fractal up to the present observational limit 00]. A crucial prob- 
lem, which is also controversial, related to the above question is whether the 
angular projection of a fractal embedded in a 3-dimensional (3d) space is 
homogeneous. Since before the extensive redshift survey the galaxy catalogs 
were for angular coordinates and now there is much angular information [Q], 
a difference in angular projection will help us choose between the two al- 
ternative pictures. Based on a numerical simulation on a fractal structure 
geneated by a Levy flight in 3d, it was stated that the angular projection 
shows power-law correlation at small angles but becomes homogeneous at 
large angles 0]. However, as we pointed out in a previous work their 
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explaination is inconsitent, taking the implicit assumption of angular homo- 
geneity in distribution. We distinguished the concepts of angular distribution 
and angular projection, the former refers to the number of points as a func- 
tion of polar angle within a given radial depth while the latter is defined 
through the solid-angular density represented as a function of polar angle. 
It was shown that the power-laws for angular projection at small angles are 
manifestinations of the so-called local angular fractal, which refers to that for 
a small conic part of a sphere defined by the angle 26 and the depth L, when 
9 is small enough, the points on each spherical shell within the conic part 
have fractal distribution. The local angular fractal may reconcil the evidence 
of fractal and that claimed to be of homogeneity, and realize Cosmological 
Principle in a fractal structure in a rather strict way. The behavior at large 
angles remained open there, and will be obtained here through the applica- 
tion of conditional calculus. Methods of calculus is also expected when one 
deals with the fractal structure on a curved space, which should be taken 
into account for large-scale structure of the universe. 

To define a function on a fractal, there are two alternative sets of indepen- 
dent variables, one is the continuous variables, or in geometric interpretation, 
coordinates of the continuous space on which the fractal is embedded, another 
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is the variables or coordinates restricted on the fractal, the inner coordinates, 
effectively they form a continuous subspace. A fractal is only a fractal looked 
in the embedding space, so we do not know whether the latter viewpoint is 
useful, practically one takes the former. 

Let us start with the derivative of /(x), f'{x) = \imAx^o[f {x + Ax) — 
f{x)]/Ax. For it to be meaningful, f{x + Ax) should be within the fractal 
set, therefore Ax should conditionally approach so that it is larger than the 
lower cut-off size of the fractal structure. Therefore we have the following 
condition of infinitesimal. 

Condition. The infinitesimal limit of change of variable should be taken in 
such a way that it is larger than the lower cut-off size of the fractal structure. 

A mathematically-ideal fractal is infinitely iterated and there is no lower 
cut-off, the above condition is automatically satisfied. Practically, this means 
that the fractal under study should be large enough. The derivative on a 
fractal has been used, e.g., in defining the conditional density [^], where the 
above condition is actually an implicit assumption. So there is almost no 
difficulty in using derivative. 

The difficulty lies in integral. How can we integrate a function /(x) 
only over the points belonging to the fractal? In the case of 1-dimensional 
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embedding space, it can be done through change of variable. Suppose Xf = 
CxX^"", where the subscript "f" denotes the variable on the fractal, is a 
coefficient and D^. is the fractal dimension, then dxf — CxD^x^'^'^dx, the 
integral is thus 



where a and b are the limits of the integration. 

Now we entend it to multiple integral, which can only be done by trans- 
formed to an iterated one. In the ordinary case of continuous space this 
is based on expressing the region of integration in terms of the infinitesi- 
mal changes of the independent variables, in geometric interpretation, the 
element of volume is expressed in the infinitesimal changes of curvilinear 
coordinates. On a fractal, the element of volume can be expressed as the 
product of the variables on the fractal, then transformed to those of the 
embedding continuous space. Here the following assumption is needed. 

Assumption. The infinitesimal element of volume on a fractal is a product 
of infinitesimal changes of lengths in orthogonal directions, these infinitesimal 
changes are fractal subsets of the infinitesimal element of volume. 

This can be called the assumption of local decomposabihty. It can be 
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seen that this assumption is based on the condition of infinitesimaL It is 
most hkely that under this condition the above assumption can be vahd, 
especially in many random fractals such as that generated by a Levy flight 
and the galaxy distribution on the scales exhibiting fractal structure. The 
sum of the dimensions of the fratal subsets is, of course, the total dimension of 
the fractal. But we have no reason that the the decomposition is independent 
of the position, the dimensions of the subfractal may depend on position. If 
they are independent of the position, the fractal can be referred to as being 
uniform. 

Thus the multiple integral on a fractal can be done as an iterated in- 
tegral over the subfractals, with the variables then changed to those of 
embedding continuous space. For example, the volume integral in terms 
of the Cartesian coordinates is / / Jyf{x)dVf — J J J f{x)dxfdyfdzf — 
CxCyCz J J J f{x)x^'^~^y^y~^z^''~^dxdydz. In this way all calculus can apply 
to fractals. 

Formally, calculus may apply with the variables of fractals by transform- 
ing derivative with the continuous variable to that with the inner variable 
of the fractal, i.e., df{xf)/dxf — c~^x^~^'^df[xf{x)]/dx. For expample, this 
derivative of number of points is a constant, showing the fractal subset form 
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an effective homogeneous space. Seen from the embedding space, this deriva- 
tive is fractionaL 

Now we turn to the angular projection of a conic part defined by L and 
29, of a fractal in 3-dimensional space. In terms of spherical coordinates 
(r, 6*, 0), the number of points in this volume is 



where A = CrCgCtf,DrDg{27r)^'t' , Dr, Dg, Dfj, and D are fractal dimensions of 
the infinitesimal subfractals and the total fractal, respectively, dr-, dg and 
dfj) represent partial derivatives. The second equality is valid when Dr, Dg 
and D^ are constants. When 6 is small, we obtain N{L,26) = AL^O^^'^^'i' , 
which is just the "local angular fractal" discussed in If an arbitrary 
azimuthal angle is considered, the number of points are then proportional 
to 0^"*. We suggest this be tested for galaxy distribution. 

For an isotropic fractal, all directions are equivalent each other, so the 
dimension of the subfractal in each direction is D/3, i.e. Dj. = Dg= D^ = 
Dx = Dy = Dz = D/3. This is most likely satisfied by the galaxy distribu- 
tion where D « 2 while « 1.3 § §. 
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The solid angle is 27r(l — cos0), so the conditional angular density defined 
from the origin is r(^) = dN {9) / dVt{e) = {dN / dO) / [dVt / d6) = {A/2tt)L^ F{e), 
where F{6) = 6^'>'^^{sin6)^'i>^^. If the sample size is characterized by 
29m and L, the average angular density over the sample is < n >m = 
N{9m)/^{9m), which is surely dependent on 9m- The angular correlation 
function is uj{9) = T{9)/ < n >0^j —1, which is dependent on 9 through 
r(^), and also on 9m- But since it breaks down spuriously at angles much 
smaller than 9m 007 the dependence of angular projection on 9 should be 
examined in T{9), thus just F{9). 

The dependence of F{9) on 9 can be seen from Figure 1, there is a power- 
law region at small angles, then there is a relatively fiat region, which is not 
so short since the plot is a log-log one, but finally it increases with 9. So the 
conclusion is that the angular projection of fractal is not homogeneous. 

For the fractal generated by a Levy fiight 0] is isotropic implied by the 
generation rule. For the galxy distribution, there is also much evidence of 
isotropy [0. So the dimension of infinitesimal subfractals are likely to be 
indeed independent on the position (cf. discussions in 0). Even if they 
are dependent on position, the qualitative nature of dependence of the con- 
ditional angular density on the angle is impossible to change. Though we 
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obtain the result disagreeing with the claim in [Q, it is not contradictory 
with what was actually reported, only a very small region was given there, 
and after a fairly flat region of ~ 0.3°, there can be observed an indication of 
increase. So the claimed homogeneity of angular projection at large angles 
is quite unsure from there. 

So the behavior of conditional angular density at sufficiently large an- 
gles can discriminate fractal model and homogeneous model for the galaxy 
distribution on the corresponding scale. Unfortunately, to our knowledge, 
there is still no such evidence. It should be noted that one should use r(^) 
instead of u!{9) since u!{9) decreases rapidly at angles much smaller than 9m- 
Among the interesting problems also is that to investigate in various frac- 
tal structures and the galaxy distribution whether the assumption of local 
decomposability is indeed valid, and the dependence or independence of the 
dimensions of local subfractals on the position. 
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Figure Caption: 

Figure 1. Log-log plot of F{9) = O^"'^ (sin 6')^'^"^ by setting Dq = 

0.6 and = 0.7 so that Dq + equals what was estimated for galaxy 
distribution p. 
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